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O 1 INTRODUCTION 

■ Al thouRh the de tection of diffuse radio emission dates back 
to lJanskvl (| 19331 ). the origin of the cosmic radio background 
, ■ (CRB) is still a mystery. The recent data obtained by the 
' Absolute Radiometer for Cosmology, Astrophysics and Dif- 
5^ , fuse Emission (arcade 2) has revived the interest in this 
Ci > question, detecting a diffuse background at frequencies be- 
tween 3 and 10 GHz that is more then 5 a above the 
COBE/FIRAS measurement of the temperature of the cos- 
mic microwave background (CMB) and well in excess of cur- 
rent estimates base d on radio source counts. More precisely 
(|Fixsen et al.ll2oTll ). the inferred value of the antenna tem- 
perature as a function of frequency can be expressed as 



ABSTRACT 

Recent measurements of the temperature of the sky in the radio band, combined with 
Hterature data, have convincingly shown the existence of a cosmic radio background 
with an amplitude of ~ 1 K at 1 GHz and a spectral energy distribution that is well 
described by a power law with index a ~ —0.6. The origin of this signal remains 
elusive, and it has been speculated that it could be dominated by the contribution of 
star-forming galaxies at high redshift if the far infrared-radio correlation q(z) evolved 
in time. We fit observational data from several different experiments by the relation 
q{z) ~ qo - /31og(l + z) with q^ = 2.783 ± 0.024 and j3 = 0.705 ± 0.081 and estimate 
the total radio emission of the whole galaxy population at any given redshift from 
the cosmic star formation rate density at that redshift. It is found that star- forming 
galaxies can only account for f^lS percent of the observed intensity of the cosmic radio 
background. 
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where Tcmb = 2.729 ± 0.004 K denotes the thermodynamic 
temperature of the CMB, Tr = 1.19 ± 0.14 K is the nor- 
malization of the radio background at //q = 1 GHz, and 
Q = —0.62 ± 0.04 is the spectral index of the CRB, con- 
sist ent with synch rotron emission from normal galaxies (see 
e.g. lCondonlll992l ). 

T he observed emiss ion is most likely of extragalactic 
origin (iKogut et al.ll201 j). and several candidates have been 
considered by ISingal et all (|2010l ). Radio source counts de- 
tected by current surveys, sensitive to flux densities above 
Si. 4 GHz 10 nJy, cannot explain more than ~ 10 per 



cent of the signal (iGervasi et al.ll2008l : iMassardi et al.ir2010l : 
IVernstrom et al]|201ll ). and low-surface brightness sources 
missed by these surveys may contribute, at most, an ad- 
ditional 15 per cent. Diffuse emission in regions far from 
galaxies is ruled out due to the overproduction of X-rays 
and 7-rays, so the only possible explanation is that the cos- 
mic radio background is dominated by fai nt sources below 
the threshold of 10 ^iJy (|Singal et al.ll2010l ). 



According to ISingal et all (|2010l '). radio supernovae 
make a negligible contribution, and radio-quiet quasars may 
be responsible for only a few per cent of the emission. Ther- 
mal bremsstrahlung from the hot gas in galaxy clusters has 
been sh own to contribute ab out 0.01 — 0.02 K at = 1 GHz 
(see e.g. IPonente et al]|201ll ). and the most reasonable can- 
didate to explain the bulk of the CRB seems to be the pop- 
ulation of ordinary star- forming galaxies at high redshift. 



Some authors (e.g. lOh I 1 19991 : ICoorav fc Furlanettd 

I2OO4I ) have tried to estimate the contribution of free-free 
emission from star-forming galaxies to the radio background 
by resorting to phenomenological prescriptions to relate halo 
mass and star formation activity at different redshifts. 
However, if the far infrared-radio correlation (FRC) ob- 
served for local galaxies holds at all redshifts, there must be 
a tight relation between the radio and infrared backgrounds 
iHaarsma fc Partridge 1 1 19981 : iDwek fc Barkeii |2002| ). From 
the measured intensity of the latter, one concludes that the 
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contribution of star-forming galaxies must be of the order of 
5 — 10 per cent. 

During the last years, the advances in infrared and sub- 
millimetric instrumentation have made it possible to investi- 
gate the evolution of the FRC over a large fracti on of the age 
of the U niverse, and several recen t studies (e.g. llvison et al.1 
l2010al lbl: iMichatowski et al.1 12OI0I ) suggest that the correla- 
tion is linear at all times, but the normalization is offset to- 
wards increasing radio loudness at high redshifts, boosting 
the expected signal from star-forming galaxies by a signifi- 
cant amount. In the present work, we make a quantitative 
estimate of the contribution of star-forming galaxies to the 
CRB. The prescription followed to assign radio luminosities 
as a function of the instantaneous star formation rate is de- 
tailed in Section (5] The evolution of the far infrared-radio 
correlation is discussed in Section |3] and the implications 
for the cosmic radio background are shown in Section [l] 
Our main conclusions are briefly summarized in Section [S] 



2 RADIO EMISSION FROM INDIVIDUAL 
GALAXIES 

In normal galaxies, ra dio emission i s always associated to 
star formation (see e.g. ICondonlll992l ) . Young, massive stars 
produce intense ultraviolet radiation that ionizes the sur- 
rounding medium, and thermal bremsstrahlung from these 
free electrons (often referred to in the radio literature as 
free-free emission) makes a significant contribution to the 
galaxy spectra in the few-GHz range. On the other hand, 
stars with M > 8 M0 explode as Type II and Type lb su- 
pernovae at the end of their life cycle. Supernova remnants 
are thought to accelerate most of the relativistic electrons 
in normal galaxies, and they constitute the main source of 
the synchrotron emission that dominates at low frequencies. 

Assuming a pure Hydrogen plasma with an electron 
temperature To ~ 10'' K, the free-free luminosity of a galaxy 
is approximately given by 



3.2 X 10-39 erg s"! Hz" 



1 ^ (gHz) (cm-3) (em's) (2) 



where u denotes the photon frequency, is a characteristic 
electron density, and 14 represents the tota l volume occupied 
by th e radio-emitting, ionized Hll regions (|Rubinlll968l : f0h1 
This volume is set by the condition that the number 
of ionizing photons Q emitted by the stars per unit time is 
equal to the recombination rate 



Q — QB Vb 



(3) 



with QB = 2.6 X lO^'^"^ cm^ s~' being appropriate for case- 
B recombination at 21 ~ 10* K. According to stellar pop- 
ulation synthesis models (e.g. iLeitherer fc HeckmanI [l995l : 
iMoUa et al.ll2009D . 



Q 



(4) 

1.5 X 1053 s-i Mq yr-i ^ ' 

where ^I* is the cur rent star formation rate (SFR) of the 
galaxy, assuming a ISalpeteJ (| 19551 ) initial mass function 
(IMF) between 0.1 and 100 Mq. In the end, the predicted 
free-free luminosity 



1.8 X 1027 erg s-i Hz-i Mq yr^i VGHz 



scales roughly proportionally with the instantaneous star 
formation rate. 

Computing the synchrotron luminosity from first phys- 
ical principles is much more involved, since it requires 
knowledge of the amount of cosmic rays injected by super- 
novae, their spectrum, and the conditions of the surrounding 
medium (most notably, its density structure and the inten - 
sity of the magnetic field). Observationally (|Condon | |l992l ). 
non-thermal synchrotron emission is about 10 times more 
luminous than the free-free continuum at u = 1 GHz, and 
its spectral index is close to ~ 0.7 for a broad range of star- 
forming galaxies. In addition, there is a tight correlation be- 
tween the synchrotron luminosity and the thermal radiation 
emitted by the dust in the infrared, which is powered by the 
stellar ultraviolet light and is thus another tracer of the star 
formation rate. The observed far infra r ed-radio correlation 
sugge sts (but see e.g. iLacki et al.ll201ol : iLacki fc ThompsonI 
|2010| . for a different point of view) that synchrotron emission 
is also proportional to the SFR, implying that 



Lb 



1.8 X 1028 erg s"i Hz^i M© yr-i VGHz 



(6) 



3 EVOLUTION OF THE FRC 

Since most of the contribution of normal galaxies to the 
cosmic radio background observed today is due to their syn- 
chrotron e mission, with thermal bremsstrahlung (see e.g. 
[oii1ll999l : ISeiffert et aDl201ll : IPonente et alll201lh provid- 
ing only a minor correction at the level of a few percent, 
equation ((6)1 has a crucial importance. In particular, the in- 
tensity of the CRB is extremely sensitive to the evolution in 
time of the relation between SFR and radio luminosity. 

It is not clear, though, whether the far infrared-radio 
correlation should evolve with redshift, and current obser- 
vational evidence i s far from being c onclusive. While sev eral 
recent studies (e.g. Ilbar et al.|[2008l : ISargent et al.ll2010l ) are 
consistent with no evolution in the FRC, some others report 
systematic trends with redshift (e.g. [V l aliakis ct al. 200 ^: 
Beswick et"al]|2008l : ISevmour et al.ll2009i : iMichalowski et al l 



2010t ). 



The main source of uncertainty is that normal galaxies 
are rather faint in the radio band. According to equation 
only the most intense starbursts, with instantaneous SFR in 
excess of ^I" > 30 Mq yr~^, would be detectable at 2 ^ 1 by 
current surveys, whose sensitivity at 1.4 GHz is of the order 
of ~ 10 /iJy. 

One pos sible solut ion (see e.g. Marsden et al.l |2009| : 
IPascale et al] [2009,; Patanchon et al.1 l2009t l is to stack the 
confusion-limited and sensitivity-limited radio images at the 
positions of thousands of infrared-selected galaxies. In do- 
ing so, one increases the signal-to-noise ratio and reduces 
the contribution of radio-loud active galactic nuclei (AGN), 
probing a population that is more representati ve of normal 
galaxie s. This procedure has been applied by llvison et al.l 
(|2010ah to a mid infrared-selected sample of galaxies, ob- 
taining that 



, iiR/3.75 X 10'" W , .-y 



Ll.4GHz/W Hz- 



(7) 



with 7 



-0.15 ± 0.03. Both the total infrared luminos- 
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4 THE COSMIC RADIO BACKGROUND 




Figure 1. Evolution of the far infrared-radio correlation. Least- 
squares fit (solid li ne) with l a limits (dash lines) to the obser- 
vationa l data from llBel]||2003l. fu ll squares), jMurphv et al.ll2009l. 
stars ), jMichalo wski et aL[[2M! p|. ope n diamonds). l lSar'gentet'^ 
20101 triangles) , llBourne et al.l2qill. full diamonds), Ijlvison et al. 
2010al . crosses) and Uvison et al.ll2010bl . open squares). 



ity LiR (defined from 8 to 1000 /imj3 A similar analy- 
sis (jivison et all l2010bl ) is consistent with no evolution, 
7 = —0.04 ± 0.03, but discarding the least reliable data at 
z < 0.5 yields 7 = -0.26 ± 0.07. 

Alternatively, one may detect high-redshift star-forming 
galaxies by observing their rest-frame infrared dust emis- 
sion, shifted towards sub-millimeter wavelengths. Based on 
a sample of 76 sub-millimeter galaxies wi th measurements in 
the radio band. lMichalowski et al.l (|2010l) conclude, in agree- 
ment with previous studies (e.g. Kovacs et al.ll200^ : lMurphvl 
I2OO9I ) that the radio emission of high-redshift galaxies scales 
linearly with the SFR, but the normalization is about a fac- 
tor of two higher than for local samples^ 

Although selection effects (see e.g. ISargent et al.|[2010l ) 
and potential bia ses arising from spectral templates 
I Bourne et aLllioill ) cannot be completely excluded, a com- 
bination of different data sets is fairly well reproduced by 



g{z) = qo~ piog{l + z) 



(8) 



with qo = 2.783±0.024 and /3 = 0.705±0.081 (see Figure^. 
Assuming that Lm oc >]/ and that the constant of propor- 
tionality does not vary with redshift, this implies that the 
synchrotron luminosity of a given galaxy scales as 



1.8 X 1028 erg s"! Hz" 



M 



yr" 



-1 (ghz) 



(l + <(9) 



In other words, we assume that the infrared luminosity is 
an unbiased tracer of the SFR and that all the evolution of 
the FRC is due to the conversion between SFR and radio 
luminosity. 



^ The difference with e.g. the far-infrared band (from 60 to 
100 ^m) is about a factor of two. and the radio power Li ,4 gHz are 
given at the rest-frame of the source, using a k-correction based 
on spectral templates. 



The specific intensity of the cosmic background at any given 
frequency is given by the integral along the line of sight 



^u'{z) 



[l + z)E{z) 



Az 



(10) 



of the average emissivity per unit volume e^i . In this formula, 
c and Hq denote the speed of fight and the Hubble constant, 
respectively. 



E{z) = ^Jn^{l + 2)3 + nk(l + zY + Qa 



(11) 



reflects the cosmological expansion, and u' — v{l -\- z) is the 
initial frequency at which the photons observed today with 
a frequency v were emitted. We adopt a WMAP7 (seven- 
year observation) cosmology with t^ m = 0.27, Q.]^ — 0, H a = 
0.73, and Ho = 71 km s"^ Mpc"^ ([Larson et al.ll201ll ) and 
compute the brightness temperature of the CRB as 
2 



T{v) = 



2ku^ 



(12) 



using the Rayleigh- Jeans approximation, where k is the 
Boltzmann constant. 

By definition, the average emissivity at a given redshift 
is the sum 



n(*,z) L^/(1') d* 



(13) 



of the contributions of all the galaxies at that redshift, with 
n{'i',z) representing the number density of galaxies with 
SFR between \i/ and ^ -I- d^' at redshift z. As long as the re- 
lation between luminosity and instantaneous star formation 
rate is linear, L„/ — k(v' , z) ^P, as indicated by equations ((5)1 
and one can express the total emissivity 

/•oo 

e^,{z) = Hi{u',z) n(*,z) * d* = k(!/',2) p.(z) (14) 



in te rms of the cosmic SFR density p* (|Dwek fc Barker! 
l2002h . The emissivity of the free-free and synchrotron com- 
ponents can be taken into account simultaneously as 



e^,{z)= [ks{u' , Z) + Ksyni)^' , Z)^ p,{z) 

with 



Kff (j/', Z) 



1.8 X 1027 erg s-iHz-iMgVr 



(15) 



(16) 



and 



/^syn 7 2) 



1.8 X 1028 erg s-iHz-iMqVi- 



{1 + zY (17) 



where [3 = Q for a non-evolving far infrared-radio correlation, 
and j3 = 0.705 ± 0.0801 to fit the data plotted in Figure [T] 
The evolution of the cosmic SFR density has been 
extensively studied during the last decade, and several 
compilations of observational d ata exist in the literature 
(e.g. [ Somcrvillc ct al. 20 01; A s casibar et al. 2002; Hopkini 
l2004l ; lHopkins fc Beacomli2006i : iMichalowski et al.ll2010l ). In 
the presen t work , we have adopted the parametrization of 
ICole et all (|200ll ) 

p*{z) 



Mq yr-1 Mpc- 



a + bz 



(18) 
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0.0137 K VGHzy ^^^^ 
whereas, using expression (|17|) . synchrotron emission yields 



0.0817 K VGHz 
for /3 = and 

^svn 



0.1402 K 



GHz 



(20) 



(21) 



for 13 = 0.705. As can be seen in Figure [S] the signal is 
dominated by galaxies at 2 < 3, due to the combined effects 
of distance dimming and the declining behavior of the cosmic 
star formation rate at high redshift. 



Figure 2. Cosmic star formation history. The solid line shows 
the best fit provided by expressio n JTSj to the data points com- 
piled bv lMichalowski et al.f ll2010lh and dotted lines illustrate an 
uncertainty of a factor of 2. 
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Figure 3. Integrated radio emission, observed at 1 GHz, from 
normal galaxies up to redshift z. AT refers to the excess signal 
above the CMB temperature. The solid line shows the contribu- 
tion of synchrotron emission, assuming /3 = 0.705, and the errors 
associated to the least-squares fit (~ 0.01 at the end of the in- 
tegration) arc indicated by the dashed lines. The contribution of 
free-free emission is plotted as a dotted line. 



and fit the observational data points reported in Table A. 4 of 
iMichalowski et al] (|2010l 'FI. The optimal values {a,b,c,d) — 
(0.011, 0.097, 2.73, 3.96) have b een obtained by means of 
the FiEstAS sampling technique (|Ascasibaij|2008l ) . a Monte 
Carlo integration scheme base d on the Field Estimator 
for A rbitrary Spaces (FiEstAS: lAscasibar fc BinnevI l2005l : 
lAscasibar.,2oTol '). The resulting cosmic star formation his- 
tory is plotted cis a solid line in Figure (2] 

Combining expressions ([TOj, gUl, and f\TE\l . 

we estimate that the contribution of free-free emission from 
star-forming galaxies to the cosmic radio background is 



^ Conversion to a Salpeter IMF between 0.1 an d IOOMq and a 
WMA P7 cosmology (following the prescription in lAscasibar et aU 
I2OO2I) amounts to a negligible correction. 



5 DISCUSSION AND CONCLUSIONS 

In the present work, we have considered the relationship 
between the cosmic star formation rate and the radio back- 
ground from star-forming galaxies in the light of recent mea- 
surements of the far infrared-radio correlation at different 
redshifts, attempting to give a further look into the signif- 
icant missing flux that has been reported by the ARCADE2 
team. 

Our main result is that normal galaxies can not be re- 
sponsible for the observed signal. 

Although we think this conclusion is fairly robust, there 
is always some room for uncertainty. Radio emission in local 
galaxies has been t horoughly stud ied, and its properties are 
well known (see e. g. I Condonll 1 9921 ) . but different gas compo- 
sitions and/or temperatures may affect the conversion factor 
between SFR and radio emission by a significant amount, of 
the order of several tens of percent. 

On the other hand, we use the cosmic star formation 
rate d ensity to constra i n the average emissivity of the uni- 
verse (|Dwek fc Barker! |2002| ). In contrast to radio source 
counts, where a population of faint objec ts below the detec- 
tion threshold is very d ifficult to rule out (|Singal et al.ll2010l : 
IVernstrom et al]|201ll h it would be extremely unlikely that 
our proposed fit underestimates the average SFR by more 
than a factor of two (dotted lines in Figure (2]). Uncertainties 
in the IMF cancel out with the production rate of ionizing 
photons given by expression Q and are not expected to 
affect the present analysis significantly. 

The most important source of uncertainty is the possi- 
ble evolution of the far infrared-radio correlation. Current 
observations seem to be compatible with /3 = 0.705 ± 0.081, 
increasing the expected emission from normal galaxies by 
about 70 per cent with respect to the case of no evolution. 
Using an extreme value P — 1 would boost the signal by 
only an additional 35 percent. 

Nonetheless, it is worth noting tha t the high-redshift 
points in Figure [1] (e.g. in the samples of IMichalowski et al.l 
l2010l : lMurphv et al.|[2009l ) are dominated by sub-millimeter 
galaxies. There is some discussion in the literature that 
these sources, whose contribution to the total SFR at z ~ 
1 — 2 is only of the ord er of ten per cent (see Figure 4 in 
IMichalowski et al.|[201ol ). may be radio-bright compared to 
normal galaxies and introduce some evolution in the ob- 
served FRC that does not apply to star formation as a 
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In fact, one would expect on theoretical grounds 
that the FRC of normal galaxies evolved in the opposite di- 
rection (/3 < 0). On the one hand, star formation at z ~ 1 
is heavily obscured by dust, and the approximation that all 
the ultraviolet luminosity is re-radiated in the infrared is 
very good. In the local universe, some fraction of the ion- 
izing photons is able to escape, and the infrared luminosity 
per unit SFR should be lower. On the other hand, galaxies 
at high redshift should produce less radio emission because 
the energy density of the CMB scales as (1-f z)**, and the rel- 
ativistic electrons injected by supernovae los e more energy 
through inverse Compton scattering (see e.g . ICarilli fc YunI 
ll999l : ICariUi et al.ll2008l : iMurphv et al.ll2009l ). Both effects, 
especially the latter, would only strengthen our conclusions, 
and the estimate with /3 = 0.705 should arguably be re- 
garded as an upper limit. 

According to our results, radio emission from star- 
forming galaxies could explain up to ^ 13 per cent of the 
intensity of the CRB. Even taking all the possible uncer- 
tainties into account, we are still far from the 1.19 Kelvin 
reported by ARCADE2 at 1 GHz. Although evolution of the 
FRC at 2 < 3 has to be further investigated, current data 
strongly suggest that it only results in a relatively minor 
boost to the contribution of normal galaxies, and hence we 
can rule them out as the main source for the radio back- 
ground. As shown in Figure |3l the contribution of galaxies 
at higher redshifts is negligible. 

Since relatively bright point sources, as well as Galactic 
or extragalactic di ffuse emission have also been ruled out 
llSingal et al.l |2010| . and references therein) , there are few 
alternatives left to explain the observed cosmic radio back- 
ground. Some possibilities are: 

(i) The ARCADE2 measurement is incorrect, or it is con- 
taminated by Galactic foregrounds. Being perfectly con- 
sistent with independent measur e ments at longer wave- 
lengths (e.g. [ Haslam et al.l Il982l : Ide Oliveira-Costa et al.l 
I2OO8I : [Rogers fc Bowrnan 20081 ). we think this possibility is 
unlikely. 

(ii) Faint star-forming galaxies at high redshift are ex- 
tremely radio bright, perhaps due to an enhanced magnetic 
field or AGN activity with respect to the brightest objects a t 
that redshift (the possibility favored bv lSingal et al.ir2010l ). 

(iii) There is a new population of numerous and faint ra- 
dio sources waiting to be discovered. 

To sum up, the nature of the cosmic radio background 
poses an exciting challenge for radio astronomy, to be faced 
in the upcoming era of Expanded Very Large Array (EVLA) 
and Square Kilometre Array (SKA). 
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